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Abstract. - We study the effects of externally applied shear flow on a model of suspensions 
of motors and filaments, via the equations of active hydrodynamics [PRL 89 (2002) 058101; 
92 (2004) 118101]. In the absence of shear, the orientationally ordered phase of both polar 
and apolar active particles is always unstable at zero-wavenumber. An imposed steady shear 
large enough to overcome the active stresses stabilises both apolar and moving polar phases. 
Our work is relevant to in vitro studies of active filaments, the reorientation of endothelial cells 
subject to shear flow and shear-induced motility of attached cells. 



The collective behaviour of suspensions of active particles( 1 ) is fundamentally different 
from that of their passive counterparts, as highlighted by recent studies on hydrodynamic 
instabilities [1,2], rheology [3] and steady state patterns [4-6] of active particle systems. Recent 
experimental studies on model systems such as motor-filament extracts [7, 8] and suspensions 
of swimming bacteria [9-11] or spermatocytes [12] have also reported strong departures from 
equilibrium behaviour. 

One surprising prediction [1] is that, unlike in passive suspensions, long-range uniaxial 
orientational order in active Stokesian( 2 ) suspensions of polar particles is always destroyed by 
a hydrodynamic instability. This is because the flow generated by a small, long-wavelength 
splay or bend deformation imposed on an oriented configuration always acts instantaneously 
(in the Stokesian regime) to further increase the deformation [1,13]. 
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( x )By active particle we mean an object which absorbs energy from its surroundings and dissipates it in 
the process of carrying out internal movements; examples include self-propelled organisms, motor-filament 
complexes, and agitated monolayers of granular particles. 

( 2 )wherc viscosity dominates and inertia is ignored 
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If only this instability could be controlled, we would have an orientationally ordered phase 
of active particles with rather unusual rheological properties [13]. We already noted two 
dramatic rheological consequences [2,3]: (i) the orientationally ordered phase is characterised 
by a nonzero steady-state average of the deviatoric stress - a kind of yield stress - and (ii) on 
approaching the orientationally ordered state from the isotropic fluid, a suspension of active 
contractile elements exhibits solid-like behaviour without translational arrest. However to 
obtain this novel material we would first need to stabilise this orientationally ordered phase. 

In this paper wc ask whether an imposed shear flow can suppress this instability of active 
Stokesian suspensions. We address this question by constructing a set of coarse-grained equa- 
tions for the hydrodynamics of the active particle suspensions. We consider apolar (nematic) 
as well as polar order, keeping in mind that the latter correspond to phases with a nonzero 
drift velocity of the active particles. Here is a summary of our results, valid in both 2 and 3 di- 
mensions, (i) The instability obtained in [1] occurs for apolar ordered suspensions as well; (ii) 
Imposed shear larger than a critical strain rate j c stabilises the orientationally ordered phase 
yielding a stability diagram (Fig. 4) controlled by two variables, a flow alignment parameter A 
and the the ratio of the shear to active stress. The drift velocity of the vector-ordered phase 
depends on 7 and A. Possible experimental tests and some complications regarding boundary 
conditions and back-flow are discussed at the end of the paper. We now turn to the details of 
our calculations. 

(A) Hydrodynamic equations : The hydrodynamic variables for a suspension of active 
particles [1,3] are the total momentum density g = pu of the particles + fluid, where u is 
the hydrodynamic velocity field and p the density, the concentration c(r, t) of active particles, 
and the orientational order parameter field. For polar systems, the order parameter is the 
polarisation vector p(r, i) of the force dipoles associated with the active particles. Active 
systems being out of thermal equilibrium, a local polarisation always implies a local average 
drift velocity of the active particle relative to the fluid. We thus take p to be the velocity field 
of the active particles relative to the fluid. Apolar nematic orientational order is characterised 
by a traceless symmetric tensor Q = S(pp — p 2 (l/d)I), where p is now the director and S is 
the scalar nematic order parameter [14]. We focus on the polar case, and indicate differences 
where appropriate for the apolar case. Momentum conservation in the Stokesian limit means 
V • a = 0, where the total stress tensor cr = cr a + cr r + <r d is the sum of contributions from 
activity, order-parameter gradients and viscous dissipation, and overall incompressibility of 
the suspension means V • u = 0, p — po, the mean density. 

The equations for the active polar order parameter p read 



In 0}, the first three terms on the left are the material derivative of p (co- moving 
and co-rotating with the suspension) and the square brackets (as well as the £ term on 
the right) contain symmetry-allowed polar contributions [15] (of which the Ai term alone 
[2] is active), ruled out in apolar nematohydrodynamics, whether passive [14, 16] or ac- 
tive [1]. The first term on the right, together with the V x u term on the left, lead to 
flow- alignment [14]. The relaxation dynamics is contained in the order-parameter molecular 
field h = c(r,t) lap — /3|p| 2 p + A"V 2 p] which favours a fixed length for p and assigns an 
clastic cost to inhomogeneities in the one-Frank-constant approximation. 

Conservation of active particles states 



d t p + (u • V) p - - (V x u) x p + [Ai (p • V) p + A 2 (V • p) p + A 3 V|p| 2 ] 




(1) 



d t c 



V- [c(u + p)]. 



(2) 
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Fig. 1 - Suspension of active vector ordered particles subject to an imposed velocity along the A-axis, 
Uo = ^yx with a gradient along the y-axis. The shear alignment angle 9 and the angle £ between the 
wavevector q and the ordering direction are also shown. 



The reactive stress from order-parameter gradients is 

<T^-^(ph + (h P ) T )+ni (3) 

where II is a generalised pressure; and the viscous stress is cr d = ^ (Vu + (Vu) T ) = rj A, 
where r\ the bare shear viscosity of the suspension. 

To determine the active stress, we make use of the fact that the simplest active particle, 
on long timescales, is a permanent force dipole [1,3]. To leading order the deviatoric part of 
the stress coming from activity is given by 

cr a (r > t) = Wc(r,t) (pp-p 2 ±l) , (4) 

where d is the spatial dimension, and the magnitude and sign of W characterise the nature of 
the elementary force dipoles [1,3]. 

Here we use this hydrodynamic description to study the effect of shear flow on the ordering 
and stability of active particle suspensions. 

(B) Shear flow : We impose a planar shear flow (Fig. 1) along the A-axis, with a velocity 
gradient along y, giving rise to an imposed velocity field uo = jyx. 

The instability of the orientationally ordered phase [1] and its stabilisation due to shear 
is most simply seen in two dimensions (d = 2); we therefore present detailed calculations in 
d = 2, and merely state results in d = 3. 

We look for homogeneous steady states of CEJ, J3J); the steady state concentration is a uni- 
form c(r, t) = cq. The amplitude and phase of the steady state active vector order parameter 
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are, respectively, 



tea, = y/^J. (6) 

The flow alignment parameter A can take values between 1 and oo, corresponding to 
< < 7r/4. Note that while p increases with shear rate, the phase 8 is independent of it. 

( C) Stability of orientational order : To determine the stability of this homogeneous steady 
state, we set p(r,i) = po(cos#, sin#) + 5p(r,i), u(r,i) = u + 5u(r,i) and c(r,i) = c + 
<5c(r,t), where the perturbations are assumed small. It is convenient to decompose <5p and <5u 
parallel and perpendicular to the ordering direction po, e.g., <5p = (5p\\, Sp±), and to eliminate 
the "massive" field 8p\\ in favour of the remaining variables yielding, after a spatial fourier 
transform, 



i 1 7 O 

d t + iq\\ Aipo + 2iq±F\ 3 p + 7V A 2 - 1 Sp± (q) + - (F cos 8 - sin 0) — (q±Sp± (q)) 



dq y 



2po9|| 



(1 + Acos2£) 



Sg±{q) + 



Sc„, 



(7) 



where P = 2Ta + 7 (A 2 — l) , F = -y/P and D = (/P. Since the instability we are interested 
in has maximum growth rate [1] at zero wavenumber, we ignore the term qi-^p operating on 
the fields, whose effect can be shown to vanish for q — > 0. 

Fluctuations in the hydrodynamic velocity field, <5u = 5g/ po, to linear order are governed 
by the Stokes equation. Imposing incompressibility and eliminating the pressure, eliminating 
8p\\ in favour of the remaining fields as above, and fourier-transforming in space, we find 
q\\ = q-x cos 8 + q y sin 8 to obtain, 







'iWqffli 



SCn 



2iWa 



m\\q± 



paq 



q±Sp±{q) _ l J^Sl S p ± (q). (8) 



Concentration fluctuations to linear order are given by 

d t Sc q = -i po5c q - ic (q ± - Fq\\ ) Sp ± (q) . 



(9) 



We are now in a position to compute the full fluctuation spectrum. While we have analysed 
the linear stability of the orientationally ordered phase in the basis spanned by Su, Sp and 
<5c, we find that the origin of instability and its recovery by the shear flow can be understood 
even in the absence of the concentration equation. Thus, to make the subsequent analysis 
more transparent, we drop the concentration fluctuation terms at the outset. 

Following [1], we express the fluctuation spectrum in terms of the splay fluctuation $ = 
V±5p±, the in-plane expansion rate 6 = V±5g±, and the angle £ made by q with the ordering 
direction: 



d t + iq\\XiPo + 7V 7 A 2 - 1 ^1 - + 0(7 



ipm 



and 



= -m& q 



2poq\\ 
iWcoq\\ cos 2^ 



(1 + A cos 2£) 



Pa 



1>„ 



O q = (10) 
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Fig. 2 - The growth rates f2±(q = 0) versus 2£ (angle between perturbing and the ordering direction), 
for either sign of activity W . Inset shows the induced flow lines arising from active rods subject to a 
flow. 



This implies that the splay fluctuations $ g have a growth rate 



First recall the generic instability when 7 = 0: the oriented phase is always unstable, either 
to splay or to bend fluctuations, depending on the sign of W [1]. This can be seen even at 
<7 = 0, where the growth rate f2+(0) > for -tt/4 > £ > n/4 when W > 0, and S!_(0) > 
for 7r/4 > £ > 37r/4, when W < (Fig. 2). The dispersion curve 0± is shown in Fig. 3 
: fluctuations with wavenumber q < qo oc \W\ grow in time; this sets the scale over which 
orientational order can be stabilised. Note that both the polar phase, as noted by [1], and the 
apolar nematic as well, are generically unstable in the absence of shear. This is clear, since 
the polar terms containing {Xi} and £ in Q, do not appear in the growth rate equation 

The imposed shear flow represented by a horizontal line 7 = const, in Fig. 3, cuts this 
dispersion curve at q* (7), suggesting that fluctuations whose scale is smaller than 1/q* are 
the first to be stabilised by the shear flow. Our estimate of q* is qualitative : the crossover 
from f2-dominated to 7-dominated at nonzero q cannot strictly be evaluated in our treatment, 
since we dropped the q x dq y terms on the grounds that they wouldn't matter at q — where 
the instability is fastest. As 7 is increased, q*(j) decreases, till at 7 = j c , this cutoff scale 
moves to zero, as q* — (7 — 7c) 1 / 2 - At this shear rate j c , the oriented phase is completely 
stabilised by the shear flow. This defines a stability boundary as a function of A and W. 

The stability phase diagram is best represented by defining a dimensionless active Peclet 
number, Pe a = 2r]j/\W\co, as the ratio of the imposed shear rate to the typical shear-rate 
produced around the active particles. As one crosses from the unstable to the stable region in 
the plane of Pe a and flow-alignment parameter A (Fig. 4) the orientational order parameter 
sets in at the value given by Eq. which in effect gives a discontinuous transition since the 




(12) 
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Fig. 3 - Dispersion curves fi± versus q, for £ corresponding to maximum growth rate. The horizontal 
line corresponding to a frequency 7VA 2 — 1 (l — ^j) defines the scale q* beyond which fluctuations 
are stabilised. As 7 increases towards j c , the stabilisation scale goes to zero as (7 — jc) 1 ^ 2 - 



order parameter in the hydrodynamically unstable region is zero. For polar active particles, 
the shear-stabilised oriented phase has a nonzero drift of the particles with respect to the 
solvent. 

Note that the critical shear rate required to stabilise the oriented phase j c is larger for 
positive W than negative. This is consistent with our earlier observation [3], that the flow 
induced by the active stresslets, oppose the imposed flow when W > and enhance it when 
W < 0. 

(D) Shear flow induced stabilisation in d = 3 : The calculation of the stability diagram is 
more tedious in d = 3; however since the spirit is the same, we merely quote results and point 
out differences. To start with, we note that po, the steady state orientational conformation 
in the presence of the shear flow Uo = ^yyx is the same as ©, i.e.. it still lies in the plane 
of the imposed velocity and its gradient (the xy plane). We then decompose the fluctuations 
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in an appropriate orthonormal basis, viz., 5p = (5py, 6p±, 5p z ), where 5pu — Sp ■ po, 5p± = 
8p- (z x po) and 5p z — Spz (similarly for <5g). As before we ignore concentration fluctuations. 

Once again, 5p\\ is massive, and we rewrite the order parameter fluctuations in terms of the 
in- plane splay <f>j^ = V ±Spj_ and <$> z = V z Sp z . Similarly invoking incompressibility, we rewrite 
the momentum fluctuations as the in-plane expansion rate 8j_ = V±5g± and 9 Z = V z 5g z . 
Eliminating the momentum fluctuations via force-balance, we find that the resulting linearised 
dynamical equations for <£>^ and $ z give rise to an eigenvalue spectrum qualitatively resem- 
bling Thus even in d = 3, one may define a stability phase diagram in the Pe a — A plane; 
a large enough shear rate stabilises an orientationally ordered phase. 

Note that the vector-ordered phase has a non-zero drift velocity Vd, with a magnitude 
proportional to cqPo (©• This macroscopic particle current will result in a counter solvent 
flow of the same magnitude. 

Experimental realisations in a planar shear flow geometry will necessarily have to contend 
with finite boundaries; it is therefore important to specify boundary conditions for the active 
order parameter p. This is especially important in the polar case, since the local polarisation 
implies a local average drift velocity of the active particle relative to the fluid. Assuming 
the active particles cannot penetrate the walls, p must be tangent to the confining walls (the 
homogeneous boundary conditions of liquid crystal physics). In plane Couette flow, however, 
this is at odds with the flow- alignment requirement of p pointing at an angle to the suspension 
velocity as in Eq. ©. This conflict must be resolved by a boundary layer in p at the walls. 
The case of wall- normal (homeotropic [14]) alignment will be discussed elsewhere. 

In conclusion we have shown how to stabilise the orientationally ordered phase of an active 
particle suspension by imposing a uniform shear flow. We have determined the nonequilibrium 
phase diagram in the plane of "active Peclet number" and flow-alignment parameter. This 
sets the stage for a study of the unusual rheological features of this shear-stabilised oriented 
phase of active matter [13]. In a forthcoming submission we will use this framework to study 
the dynamics of reorientation of endothelial cells subject to shear flow [17]. 

We thank the Indo-French Centre for a grant (IFCPAR Grant # 3504-2). 
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